Collective modes of a helical liquid 
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We study low energy collective modes and transport properties of the "helical metal" on the surface 
of a topological insulator. At low energies, electrical transport and spin dynamics at the surface 
are exactly related by an operator identity equating the electric current to the in-plane components 
of the spin degrees of freedom. From this relation it follows that an undamped spin wave always 
accompanies the sound mode in the helical metal - thus it is possible to 'hear' the sound of spins. 
In the presence of long range Coulomb interactions, the surface plasmon mode is also coupled to 
the spin wave, giving rise to a hybridized "spin-plasmon" mode. We make quantitative predictions 
for the spin-plasmon in Bi2Se3, and discuss its detection in a spin-grating experiment. 
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Introduction - Recently, topological insulators have 
been theoretically predicted and experimentally observed 
in both quasi-two dimensional f2D) and three dimen- 
sional (3D) systems[l|, i, d, i, S E Q- The concept 
of a topological insulator can be defined within the non- 
interacting topological band theoryj^, Hor more gener- 
ally within the topological field theory [10|, which is also 
valid for interacting systems. The simplest topological 
insulators such as Bi2Se3 and Bi2Te3 have a full bulk 
insulating gap and a surface state consisting of a sin- 
gle Dirac cone0, H, 0|- As is the case for the helical 



edge states of a 2D topological insulatorfllL Il2l|. the spin 
and the momentum are intimately locked in the "heli- 
cal metal" surface state of the 3D topological insulator. 
This locking effect has been theoretically predicted [yl for 
Bi2Se3 and Bi2Te3, and experimentally observed[l9j in 
Bi2Se3. 

In this paper, we study the universal surface state 
properties of the simplest 3D topological insulators, and 
consider a system governed by a single isotropic Dirac 
cone at energy scales much lower than the bulk insulat- 
ing gap. We study the consequences of the helical nature 
of the metallic states: the coupling between spin and 
charge excitations, and collective modes of the helical 
liquid. Our theory here is directly applicable to the case 
of the Bi2Se3/Bi2Te3 family, which has a single isotropic 
Dirac cone that remains isotropic for low dopant concen- 
trations. 

Spin dynamics and electrical transport - Starting from 
a low energy effective Hamiltonian for the bulk of a 3D 
topological insulator in the Bi2Se3 family, the low energy 
surface Hamiltonian was derived by Zhang et al. Q by 
diagonalizing the bulk effective Hamiltonian with open 
boundary conditions, and by integratingout the high en- 
ergy bulk degrees of freedom. In Ref . [6| , it was shown 
that for a surface in the xy-plane, the helical states at 
low energies are governed by 

H ^ j d^x?At(x) [hvf{i X {-iV)) ■ cr - n] V'(x), (1) 

where k is the Bloch vector in the 2D surface Brillouin 



zone and cr the Pauli matrices describing electron spin. 
Here, /i is the value of the chemical potential relative to 
the surface Dirac point. This Hamiltonian adequately 
captures the dynamics of the surface for energies much 
smaller than the bulk gap A, and for length scales much 
larger than hvf / A, the penetration depth of the gapless 
surface states into the bulk. In this regime, there are 
several universal features that characterize the surface 
states. First there is the operator identity relating the 
charge current to the in-plane component of the spin on 
the surface: 



j(x) = ^''(x)vf (cr X z)-0(x) = u/S(x) 



(2) 



Such a simple relation between charge current density 
and spin density is the key observation of this work, 
which is a unique property of the helical liquid, and 
leads to many intrinsic correspondences between spin and 
charge dynamics in the system. It should be noticed that 
the operator identity Eq. ([2]) remains valid even when in- 
teraction terms are added to the Hamiltonian. The only 
deviation from this identity comes from the change of 
Fermi surface shape due to rotational symmetry break- 
ing, which has been observed in Bi2Te3 when the chem- 
ical potential is close to the bottom of the bulk conduc- 
tion band 7] . In this regime, the Fermi velocity obtains 
an 0{k/kfY distortion[l3|,|lJ], and the operator identity 
above no longer rigorously holds. However, there is still a 
one-to-one correspondence between the velocity and spin 
of the surface electrons, so that many conclusions we will 
discuss in the following will still hold qualitatively in this 
regime. 

In what follows, we will assume a perfectly circular 
fermi surface in the continuum limit. From the above 
identity dH), we derive the important dynamical identity 
between correlation functions 



v'jeikejl{TrSi{T)s-j{Q)) = {Trjk{T)jl{0)) 



(3) 



which relates the electric transport to the dynamical spin 
structure factor in the linear response regime. The elec- 
tric response of the helical metal can thus be calculated 
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from the generalized spin and density susceptibility 



susceptibility tensor in the basis (n, s-^) has the form 



(4) 

where tr^ = [l^cr), and we have introduced the Matsub- 
ara Green function Q. In particular, the optical conduc- 
tivity cra;a;(w) is simply related to dynamical spin suscepti- 
bility: Gxxi^) — z^'^'fXyyi^)! while the Hall conductivity 
is related to the off-diagonal dynamical spin susceptibil- 
ity a^yiuj) = —■^v'jxyxi^)- In principle, this relation can 
be verified by comparing the optical conductivity with 
the experimental results of spin-polarized neutron or po- 
larized light scattering. 

Due to the rotational invariance of the Fermi surface 
about the z axis, the longitudinal (s^ = q • s) and trans- 
verse (s"^ = z • (q X s)) components of the spin excita- 
tions are decoupled. In particular, it can be shown that 
susceptibilities involving correlations between longitudi- 
nal and transverse spin degrees of freedom vanish iden- 
tically. An explicit evaluation of Eq.(|3]) shows that the 
4x4 susceptibility tensor decomposes into two 2x2 ma- 
trices, one with the density and transverse spin, and the 
other with the perpendicular (s^) and the longitudinal 
spin. The density-transverse spin block of the suscepti- 
bility tensor can be determined from the Ward identity 
that has its origin in the continuity equation for the den- 
sity: dtfiq — —iqjq, where = q • jq is the longitudinal 
current. The Ward identity, combined with the operator 
identity in Eq.(l2]) gives us 



dtria 



■ T 



(5) 



We note that the longitudinal current is simply the trans- 
verse spin degree of freedom. 

The structure of the density and spin response bears 
a remarkable resemblance to source-free Maxwell elec- 
trodynamics in 2 + 1 dimensions, as first pointed out in 
Ref. [l5| . The analogy is made precise when we identify 
the density with the "magnetic field" perpendicular to 
the surface, and the transverse spin components with the 
"electric field" in the plane of the surface. Eq. ([5]) can be 
written a,s dtn — — V x \n the real space. The equation 
of continuity that connects the density to the transverse 
spin is thus precisely the Faraday law for this system. 
Moreover, Gauss' law for the electric field V • = is 
satisfied by construction, since the spins are transverse 
to the in-plane momentum. A similar identification has 
been made for a system with Rashba spin-orbit coupling 
[l5| : however, the analogy to electrodynamics is less pre- 
cise in Rashba systems than in the case of the helical 
metal due to the presence of the quadratic dispersion in 
the Rashba Hamiltonian. 

Applying the Ward identity to the response functions, 
we find that the density and transverse spin subset of the 
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(6) 



We note that in order to obtain this form, we need to reg- 
ularize the susceptibility tensor, taking the fi = ground 
state as the vacuum 16l. |17|. 

The bare susceptibilities can be evaluated explicitly 
for arbitrary lu and q using the non-interacting Green 
function 



^(°)(k,ic^„) = K,-7^(k)]-i= ^ 7',Gi°Hic^„,k), 
1, 



s=±l 



Vs 



[1 + ,s(z X k) • cr] 



G(") = [iUn- SVfk + fl]-^ 



(7) 



As shown in Fig. [U we find that particle-hole excitations 
are found everywhere except for Vfq < uj < 2/1 — Vfq. In 
this region, the imaginary components of all the suscep- 
tibilities vanish and collective excitations are undamped. 




FIG. 1: Dispersion of the collective modes in the long wave- 
length regime. The plasmon (solid red line) is shown for 
a = 0.6 and the zero-sound (dashed blue line) mode is shown 
for U = 0.3i7c. The shaded region denotes the particle-hole 
continuum where the imaginary parts of the density and spin 
susceptibilities are non-zero. The collective modes are long- 
lived in the unshaded region v/q < lj < 2^ — v/q. 

Collective modes - In the long- wavelength, low fre- 
quency limit, q <ti fc/, w <C fJ., the non-interacting density 
susceptibility to leading order is 



- ^ 
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a; + Z(5 
Vfq 



(8) 



with corrections 0{q/kf)'^. Thus, for uj > Vfq, the imag- 
inary part of the density and spin susceptibilities vanish 
and undamped collective modes exist. To determine the 
collective mode spectrum in this regime, we treat the 
susceptibility tensor within the random phase approxi- 
mation (RPA). 
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In the case of the long-ranged Coulomb interaction, 
U{q), the RPA correction to the susceptibilities in Eq.® 
is given by 



X'"'''(q,u.' + a) = i 



1 -1 



X 



(9) 



where x is the density and transverse spin susceptibility 
in Eq.([S]) and is a Pauli matrix in the same basis. 
Although the susceptibility is a 2 x 2 matrix, the result 
of the RPA correction turns out to be the same as an 
ordinary 2D Fermi liquid: 



For short-ranged Hubbard- like interactions, the RPA 
susceptibilities are 



X 



'RPA 



(13) 



The reason for this difference in the RPA susceptibilities 
is that short ranged Hubbard-like interactions can be de- 
composed in both the density and spin channels, both of 
which contribute to RPA diagrams, whereas the Coulomb 
interaction can only be decomposed in the charge chan- 
nel in the RPA approximation. In both cases, the RPA 
susceptibilities satisfy the Ward identity. 



, 1 X \ Xnn 



VfQ 



(10) 



Collective mode excitation spectra are obtained via the 
poles of the matrix of RPA susceptibilities. For the 
Coulomb interaction U{q) = 2Trhavf/q, where a — 
e^/edhvf {id being the dielectric constant of the topo- 
logical insulator) is the fine structure constant for the 
helical metal. Therefore the plasmon dispersion satisfies 



hv 



Vx^^ 



- 1 



(11) 



which for small q establishes that plasmon modes are 
gapless and propagate with a dispersion 




a q 
2^ 



f 



(12) 



similar to an ordinary Fermi liquid in 2D. Note that in or- 
der for the plasmon to propagate, the solution to Eq. lTTI) 
needs to stay in the unshaded region of Fig. [1] which 
means that for frequency above ujc = (ck/2)(M/^): the 
plasmon will merge into the particle-hole continuum. We 
expect Huic ~2.2meV for the 612863 family; this comes 
from ed « 100 [l3] and Vf ~ 5.0 x lO^m/s with /i typi- 
cally around lOOmeV 

In sharp contrast to the conventional Fermi liquid, the 
surface plasmon of the helical liquid always carries spin, 
and could be appropriately called the spin-plasmon. This 
is a consequence of the equation of continuity for the 
density and the operator identity, Eq. O dtn^ = —iqs^. 
Thus, a fluctuation in the density will be accompanied by 
a transverse spin fluctuation due to the "Faraday law" 
that couples them. In Fig. [2l the nature of the spin- 
plasmon is shown pictorially at a fixed oj and q. A den- 
sity oscillation induces a transverse spin wave in perfect 
analogy with Maxwell electrodynamics in 2 -t- 1 dimen- 
sions. At the peaks and troughs of the density wave, 
the transverse spin components are zero, whereas regions 
where the density variations vanish are accompanied by a 
maximum spin polarization. The locking of the spin and 
charge degrees of freedom in the spin-plasmon collective 
mode is unique to the helical liquid and marks a striking 
difference relative to the ordinary Fermi liquid. 



FIG. 2: The spin-plasmon collective mode: a density fluc- 
tuation (green surface) induces a transverse spin fluctuation 
(red arrows) or vice-versa. To detect the spin-plasmon, a spin- 
density wave is generated by a spin grating. Two orthogonally 
polarized non-coUinear incident beams (relative angle of 29) 
induce a spin polarization wave. The transverse component 
of the photon helicity is P cos 9 cos ip where P is the total pho- 
ton helicity amplitude and is the beam plane makes with 
the surface. The induced plasmon charge oscillation can be 
detected by conventional means. 

The zero sound spectra for short range interaction can 
be determined by the pole of Eq. (jl3p . which is given by 



27rft2„2 
= (x 



1 



Vx^ 



1 



- 1 



(14) 



In this case, there is an essential difference from the ordi- 
nary fermi liquid. The right hand side of Eq. (fTil) takes 
its value in the range +cxd), so that the equation has a 
solution only if A-Kfi^vj/Ufj. > 1, or U < Uc = 'i'Kh'^v'j / fx. 
The sound wave dispersion for the two limiting cases 
U '^Uc and U < Uc are 



(15) 



Physically, the disappearance of the sound mode for 
strong interaction U > Uc can be understood as a conse- 
quence of the attractive interaction in the spin channel: 
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note that the interaction vertex in Eq. ([13]) for the den- 
sity correlation Xnn has the opposite sign of the vertex 
for the transverse spin correlation xtt- In an ordinary 
Fermi hquid with repulsive short range interaction, the 
same argument leads to the damping of the spin wave, 
but does not affect the zero sound since the spin and 
charge response are decoupled (or only weakly coupled in 
systems with weak to moderate spin-orbit coupling). By 
contrast, in the helical spin liquid, the spin and charge 
responses are intrinsically coupled due to the operator 
identity Eq. so that the sound wave can be damped 
even for a repulsive interaction. 

Spin grating measurements - Our basic method of de- 
tecting the spin-plasmon mode is to excite the spin degree 
of freedom and to detect the propagating density wave 
coupled to the spin polarization wave. In other words, 
we generate transverse spin polarization and detect the 
induced density wave, which we can measure through 
spatial modulation of reflectivity. 

For generating the transverse spin wave, we propose 
a method similar to the transient spin grating (TSG) 
used in Refs. and 21. For our version of TSG, 
shown in Fig. [21 we need to have two non-coUinear 
femtosecond laser beam pulses incident on the topo- 
logical insulator surface, with z set as the surface nor- 
mal; the wave vector directions of these two beams are 
(± sin 9, — cos 6 cos (p, — cos 9 sin ip) . By linearly polariz- 
ing the two beams in orthogonal directions, we can gen- 
erate by interference alternating photon helicity in the 
direction (sin 9, — cos 9 cos (p, — cos 9 sin ip) with the grat- 
ing vector q along x; the q can be varied by changing 
the relative phase between the two interfering beams. 
For the persistent spin helix experiment, Ti:sapphire 
lasers (wavelength 650~1100nm) were used to obtain 
q = 0.34 - 2.5 x 10\m~-^ [IH; this q value would suit 
our purpose (a typical kf with /i lying in the bulk gap 
can be taken as 60 x lO'^cm"^ for Bi2Se3 or Bi2Te3 @,0])- 
So long as we keep the two beams' polarization orthog- 
onal and intensity equal, this photon helicity wave does 
not directly generate any density wave. However, the 
spin polarization wave it generates will have a nonzero 
s'^ (X cos9 cos ip for cp ^ 7r/2. After the laser pulses are 
applied, the transverse spin components s'^ propagate 
through the spin-plasmon collective mode. On the other 
hand, other spin components are Landau damped. Note 
that intensity / of the induced density-spin wave has ip 
dependence, which vanishes for = 7r/2 since in that 
case s'^ — 0. Observation of such dependence would be 
a strong evidence confirming our helical liquid theory. 

In conclusion we presented a general theory of the col- 
lective modes of the helical liquid. We derived a general 
relation between the charge current and the spin which 
is valid for arbitrarily interacting systems. We show that 
the the spin-plasmon mode propagates on the surface of 
a topological insulator, and propose a experimental set- 
ting to detect this mode. The spin-plasmon mode uni- 



fies spintronics and plasmonics, two frontier branches of 
current research, and can be used for spin transport in 
spintronics devices. 
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